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ABSTRACT 

We report the discovery of a planet transiting the star WASP-80 (1SWASP J201240.26-020838.2; 2MASS J20124017-0208391; 
TYC 5165-481-1; BPM 80815; V=11.9, K=8.4). Our analysis shows this is a 0.55 + 0.04 M jup , 0.95 + 0.03 R jup gas giant on a circular 
3.07 day orbit around a star with a spectral type between K7V and MOV. This system produces one of the largest transit depths so far 
reported, making it a worthwhile target for transmission spectroscopy. We find a large discrepancy between the v sin i+ inferred from 
stellar line broadening and the observed amplitude of the Rossiter-McLaughlin effect. This can be understood either by an orbital 
plane nearly perpendicular to the stellar spin or by an additional, unaccounted for source of broadening. 
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Numerous planets have been found since Mayor & Queloz 
(1995) and, like 51 Pegasi b, most orbit stars whose spectral type, 
mass, or size are similar to the Sun's. This occurs even though 
a few surveys have concentrated their efforts on other spectral 
types notably towards M dwarfs, such as the HARPS M-dwarf 
survey (|Bonfils et"aL||2013| ) or the M-Earth project (Nutzman 
& Charbonneau 2008). In their rarity, those planets nevertheless 
help us better understand the processes leading to planet forma- 
tion. 

The WASP (Wide Angle Search for Plan ets) survey aims to 
find transiting planets (Polla cco et al.||2006| ), and has now sur- 
veyed most of the night sky in both hemispheres. With some 
70 planets now publicly announced, this is the most efficient 
ground-based planet discovery project. Thanks to its observation 
of now more than 30 million stars of magnitude between 8.5 and 
13.5, it can pickup those rare planets that have avoided detec- 
tion by the radial- velocity surveys or even by the space-missions 
Kepler and CoRoT which have surveyed only 150 000 stars each. 
Amongst those ra re planets found by WASP is t he first gas giant 
around a 8 Scuti (Collie r~Cameron et al.|2010) and the popula- 
tion of very short period gas giants, such as WASP- 12, 18, 19 
and 43 flHebb et al.|2009l|20T0i [Hellier et al.|2009l|20TT) . 

Despite their numbers and the facility of discovering them 
(radial velocities or transits), the occurrence rate of hot Jupiters 
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* using WASP-South photometric observations, from Sutherland 
(South Africa), confirmed with the 60 cm TRAPPIST robotic telescope, 
EulerCam, and the CORALIE spectrograph on the Swiss 1.2 m Euler 
Telescope, and HARPS on the ESO 3.6m (Prog ID 089.C-0151), all 
three located at La Silla Observatory, Chile. The data is publicly avail- 
able at the CDS Strasbourg and on demand to the main author. 



orbiting solar-type stars is low. It has been estimated to be as 
high as 1.5+0.6% by Cumming et al. ( 2008) from radial v elocity 
surveys, and as low as 0.5 + 0.1% by [Howard et aL](|2012|) from 



the Kepler results, pohnson et aL] ( [2010| ) have made a case that, 
because no hot Jupiter was known to orbit an M dwarf, their oc- 
currence must therefore be lower. Not long afterwards Johnson 
et al. ( |2012| ) announced the discovery of a transiting gas giant 



around a star observed by Kepler, KOI-254, describing it as a 
'lone example [...] for some time to come , \ Approximately 300 
M dwarf systems have been searched for planets between the 
main radial velocity teams (Johns on et al.|[20T0| ). The M-Earth 
project is targeting about 3 000 (with a geometrical detection of 
only 5-10%). If the rate of hot Jupiters is but a half to a third that 
of solar type stars, there is a significant cha nce that such planet s 
have avoided detection, a point made by Bonf ils et al.| ( |2013| ). 
Knowing this rate is important since gas giant formation is per- 
ceived as less efficient because protoplanetary disc ma sses scale 
with their primary's mass as dynamical timescales do (Laughlin 



let al.|2004[|lda & Lin|2005l | Alibert et al.|2011t|Mordasini et al. 
|2012| . 

Within this context, we announce the discovery of a gas gi- 
ant transiting a late K-early M dwarf. We first describe our data 
collection, then its analysis, and finally the results we obtain. 

1 . Observations 

WASP-80 (1SWASP J201240.26-020838.2; 2MASS J20124017 
-0208391; TYC 5165-481-1; BPM 80815) was observed 5 782 
times during one season at the WASP-South facility in 
Sutherland (South Africa), by a single camera, between 2010 
May 15 and 2010 September 2 6 (Fig.[2]3). The Hunter algorithm 
( Collier Cameron et al.| [2007) found a period at 3.07 days and 
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Fig. 1. Co-added HARPS spectra on WASP-80 compared with ELODIE spectra of 61 Cyg A & B and of BD +47 44A over the Na i 
Doublet and H a . 



uncovered a transit-like signal from five partial events, with a 
large depth that nevertheless corresponds to a planet-sized ob- 
ject once the colours of the star indicated it was a potential M0 
dwarf. 

The star was catalogued for spectroscopic follow-up on 201 1 
May 09, and the first radial velocity measurements were ob- 
tained with the CORALIE spectrograph in July 2011. Thirty- 
seven spectra have been collected between 2011 July 21 and 
2012 September 12, including ten measurements obtained dur- 
ing the Rossiter-McLaughlin effect on 2012 June 19. We also 
used HARPS and acquired sixteen spectra during the transit of 
2012 September 10. Atmospheric conditions were poor with see- 
ing > 2" at the beginning of the sequence. Eight measurements 
were obtained in the nights leading to and following the transit, 
with one point badly affected by weather and excluded from the 
analysis. Radial velocities were extracted using a K5 correlation 
mask and those data also show a 3.07 day variation, in phase with 
the photometry. No such variation can be observed in the span 
of the bisector slope, or into the width of the line (Fig. [2^). This 
indicates a movement of the spectrum with time, as expected for 
an orbiting planet. 

To complete the confirmation of the system and obtain pre- 
cise physical parameters, three higher precision lightcurves were 
observed, two by TRAPPIST (Jeh in etal.|2011| ) on 2012 May 07 
and 2012 September 10 in the z-band, and one by EulerCam us- 
ing an r'-Gunn filter on 2012 July 26 (Fig. |2j:). Our data were re- 
duced and prepared for analysis in the same manner as in previ- 
ous WASP discoveries. Useful references can be found in Wilson 
|eTaLl(|2008]),|Anderson et aTJ ( |MlTT ), |Gffion et al.| ( [2009] [20T2] > 9 
and |Lendl et al.| ( |2012| ).~ 



2. Spectral analysis 

The HARPS spectra were co-added, leading to a single spec- 
trum of signal-to-noise of 75:1. Its analysis was con ducted fol- 
lowing the methods described in |Doyle et al.| ( |2013] ); the results 
are displayed in Table [T] The H a line being weak, it indicates a 
low effective temperature around 4000 K and a K7 spectral type. 
The TiO bands are also weak, typical of a late K-early M dwarf. 
The Na i D are very strong, implying a surface gravity (log g) 
around 4.6, closer to an M spectral type than a K7 (Fig.[T]). The 
equivalent widths of several clean and unblended Fe i lines were 
measured, in order to determine the stellar metallicity, evaluate 
microturbulence and confirm the estimated from H a . They 



were also fitted to estimate the broadening caused by the pro- 
jected stellar rotation velocity (v sin Macroturbulence was 
assumed to be zero since its effect is expected to be lower than 
thermal broadening (Gray 2008). We found v sin /* = 3.55 ± 0.33 
kms -1 . Because this seemed unusually large for this spectral 
type, we looked at all the stars with a similar B-V present in the 
HARPS archive: WASP-80 has the widest lines in the sample. 
There is no significant detection of lithium in the spectrum, and 
we can place an equivalent- width upper limit of 30mA, meaning 
log A(Li) < 0.0 ± 0.2. For ear ly M-type stars, lithium ca n be this 
depleted in less than 100 My ([Sestito & Randich|2005]). 



In their survey of high proper motion stars, |Stephenson 
(1986) listed the star as spectral type K5. This is, however, 
inconsistent with the results of the spectral analysis. An in- 
dependent check of the stellar temperature can be obtained 
from the InfraRed Flux Method (IRFM) (Blackwell & Shallis 
|1977| ), which has been used to determine T e s and stellar an- 
gular diameter (6) by estimating the total observed bolomet- 
ric flux from broad-band photometry from NOMAD, TASS, 
CMC14, and 2MASS. This gives T eff = 4020 ± 130 K and 
= 0.1 13 ± 0.008 mas, consistent with the T g q from the spectral 
analysis. There is no sign of any interstellar Na D lines in the 
spectra, so reddening is expected to be negligible. 



3. Results 

We applied the same fitting MCMC algorithm as described in 
|Triaud et al. ( |2011| ). The stellar mass was constrained by a 
Gaussian prior with mean and standard deviation corresponding 
to th e value obtained us ing the empirical mass-radius relation 
from |Torres et al.| ( [2010| ). The fit of t he model over the data in- 
forms us of the mean stellar density ( [Seager & Mallen-Ornelas 
2003 ) found to be 3.12 ± 0.02 p Q . When combined with a mass 
of 0.58 ± 0.05M o it gives a radius of 0.57 ± 0.02 R Q (entirely 
compatible with the Torres relation) and a log = 4.69 ± 0.02. 
Those values are also co mpatible with theore tical mass-radius 
relationships presented in |Baraffe et al.| ( [T998] ) (see Fig. [3]). 



The fit to the radial- velocities gives a reduced xi - 0.75 ± 
0.17, suggesting we may have overestimated our error bars on 
the measurements. The largest contribution to the^ 2 comes from 
the second point in the HARPS Rossiter-McLaughlin series, 
likely affected by high seeing. No eccentricity can be detected; 
we place a 95% confident upper value at e < 0.07. Similarly, if 
there is any other perturber in the system, it adds an accelera- 
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Fig. 2. a) Top to bottom: CORALIE (upright red triangles) and HARPS (inverted black triangles) radial velocities on WASP- 80 
plotted with a circular Keplerian model and residuals; below: change in the span of the bisector slope and change in the FWHM of 
the CCF. b) Phase-folded WASP V+R photometry with model and its residuals, c) Top to bottom: the two TRAPPIST z band and 
the EulerCam r'-Gunn transit lightcurves with models over plotted. The residuals are displayed in the same order below, d) Zoom 
on the Rossiter-McLaughlin effect showing CORALIE and HARPS radial velocities with the most likely model and the residuals 
from the fit. 



tion lower than 24 ms _1 yr _1 . The results of the fit are located 
in Table [T] The mass function f(m) and the log g p are directly 
obtained from fitting the data; they indicate we have discovered 
a new transiting planet. Using the stellar mass we obtain a mass 
and radius for our object and find 0.55+0.04 M jup and 0.95+0.03 
Rjup- 

The Rossiter-McLaughlin effect is marginally detected even 
though a semi amplitude between 60 and 70 ms" 1 had been ex- 
pected. This can be explained either by a highly inclined planet 
(the impact para meter b < 0.1, even when adjusting only the 
photometry (see |Triaud et al.| ( |2011| ) for details), or by an ad- 
ditional unaccounted-for broadening of the spectral lines. This 
could occur in the presence of magnetic fields producing a par- 
tially resolved Zeeman line splitting. The presence of strong 
magnetic fields ha s been r eported for a num ber of M dwarfs by 
|Morin et aL] ( |2010| ) and ponati et al.| ( [2008] ). We examined our 
spectra and found no signs of Zeeman broadening. 

Two different fits were attempted. For the first one, we as- 
sumed p = 0° + 20 and found V sin/* = 0.91 + 0.25 kms" 1 , in 
strong disagreement with the value inferred in section [2] For the 
second attempt, we chose to impose a prior on V sin /* = v sin /*, 
and as expected we find two well separated and symmetrical so- 
lutions for the spin-orbit angle: J3 = +75° + 4. We caution here 
that this angle is entirely dependent on the value of v sin /* . 

The rotation rate (P = 8.5 + 0.8 d) implied by the v sin /^ 
gives a gyrochronological age of ~ IOO^q My using the |Barnes| 



( |2007| ) relation. The presence of Ca H+K emission indicates that 
WASP- 80 may be a young active star. Despite this, we do not de- 
tect any rotational variability > 1 mmag in the WASP lightcurve 
or in a multi-epoch TRAPPIST campaign. None of the transits 
show signs of stellar spot crossings. 

Further more, using the value s of proper motion reported in 
NOMAD ( |Zacharias et al.|2004| ) (-100 + 3, -60 + 8 mas yr" 1 ) 
and the systemic velocity we observed (10.2 km s" 1 ), we com- 
puted its galactic dynamical velocities (U, V, W) = (30.0, -11.7, 
12.8) km s" 1 . Those values are well away from any of the known 
young moving groups reported in |Zuckerman & Song ( 2004 ). 
This means the gyrochronological age is not reliable. 

4. Conclusions 



Our observations and their analysis allow us to conclude there is 
an unseen transiting companion orbiting WASP-80 whose mass 
and radius are planetary. 

Even though only a few cool stars have been observed by 
WASP, our planet confirmation rate is similar to sun-like stars. 
We only followed-up 26 stars for which the Hunter algorithm re- 
turned a signal and whose colours indicate they are of K5 or later 
type (only three were classified as M0 including WASP-80). Our 
data show we have 17 blends (for eight of those the original tar- 
get was identified as a red giant), three spectrally resolved eclips- 
ing binaries and one unresolved, two potential triples (blends but 
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Table 1. Parameters of the WASP-80 system. 



1SWASP J201240.26-020838.2 
2MASS J20124017-0208391 
TYC 5165-481-1 
BPM 80815 



Filter 


Magnitude 


Filter 


Magnitude 




2MASS 1 




NOMAD 1 






J 


9.218 + 0.023 


B 


12.810 




H 


8.513 + 0.026 


V 


11.870 




K 


8.351 +0.022 


R 


11.110 




TASS4 3 




CMC14 4 






V 


11.881 +0.228 


r' 


11.358 ±0.011 


o 


I 


10.279 + 0.105 








Parameter 


Value & lcr error 


Parameter 


Value & lcr error 




from spectral 


line analysis 








Spectral type 


K7V 


Distance 


60 + 20 pc 






4145 + 100 K 


[Fe/H] 


-0.14 + 0.16 






4.6 + 0.2 (cgs) 


log A (Li) 


< 0.0 + 0.2 


0.50 


v sin i+ 


3.55 + 0.33 km s" 1 


ft 


0.3 + 0.3 km s" 1 


Mass 


0.58 + 0.05 M 


Radius 


0.63 + 0.15 R Q 





jump parameters for the MCMC 
Period 
Depth 

V V sin /* cos p 
VVsin/* sin/? 



3.0678504 j^j d 
0.02933 <™> 

n (+0.12) 
°' 48 (-0.13) 

7 o (+0.12) 
±l./o ( _o.09) 



0.17126 ^ 
0.07699 ^ 

, (+12) 



derived parameters from the MCMC 
Rp/R* 
R*/a 
log g* 
P* 

M* (prior) 

R* 
V sin 



4.689 
3 117^ 

(-20) 
<- ? (+0.05) 
U ' J 7 (-0.05) 

571 (+16 ^ 

7 1 (-16) 

3 46 (+a34i 



0.0346 
<0.07 



(-li) 



(cgs) 

Po 
M Q 

Ro 

kms~ 
' N AU 



T (BJD) 

Width 

b 

K 



f(m) 

Rp/a 

log £ P 

Pp 
M n 



2456125.417512 
0.08800 £J*> d 
0.019 p» R 
110.9 gms" 




0.425 10" 



-9 (+43) 
(-31) 

0.013183 (+39) 



M Q 



Rn 



Irl 



(-35) 

3 ^gs) 

(-39) PJ U P 
(+30) M . 
(-39) iVA J u P 

952 (+26) R 

U ' 7JZ ' (-27) "^JUP 

±75[^deg 
89.92^ deg 



3.178 
0.554 
0.554 



(-52) 

Fig. 3. Posterior probability density, output of the MCMC, show- 
ing the mass of the primary (drawn from a prior), and its radius 
(inferred from the transit signal). 1, 2 and 3 <x confidence regions 
are drawn. The over-plotted g rey line is the theoretical 0.5 Gyr 
mass-radius relationship from Baraffe et al.| ( |1998| ), in black is 
^he 8 Gyr relation (for solar metallicity). Side histograms show 
the marginalised parameters. 



' (-0.12) 

< 24 ms" 1 



yr 



Note: Mass and radius estimated using the Torres et al. ( 2010]) calibra- 
tion. Spectral type estimated from r eff using Table B.l in |Gray| ( [2008| ). 
Distance estimated using the IRFM angular diameter and the stellar 
radius. Units based on the equatorial solar an d jovian ra dii a nd masses 
taken from A llen 's Astrophysical Quanti ties ( Cox|2000]>. 1 - Skrutskie 
|et~aT] ( [20061 2 - |Zacharias eTaL] ( [2004} 3 - [Proege et ah] ( |2006| > 4 
ViZier I/304/out 



WASP-80b is a warm Jupiter when we consider its tempera- 
ture, and yet it belongs to the hot Jupiter population. Because of 
the high density of the host, and the low density of the planet, it is 
located about 3 Roche radii away from the star, just as would be 
expected if it had circularised from an earlier, more eccentric or- 
"bitlMatsumura et al. 2010 ). Although the Rossiter-McLaughlin 
effect is observed symmetrical, the planet's orbital spin could be 
severely inclined if the host star rotates as quickly as the spectral 
line broadening indicates. If no additional spectral line broaden- 
ing mechanism is discovered, then WASP- 80b would become a 
rare example of a severe ly inclined planet whose host star's T e s 
is cooler than 62 50 K ( [Brown et al.|2012t|Albrecht et al.|20l2| 
|Winn et al.|20T0| ). 



gravitationally bound), one false alarm, and tw o planets (WASP- 
43b ( [Hellier et al.||2(jTTl |Gillon et al.|[2012l ), and WASP-80b). 
A total of two planets out of 26 candidates is remarkably close 
to our mean discovery rate of 8.8 ± 1.2% planet per candidate 
( |Triaud|2011| ). In addition to the observed rarity of hot Jupiters 
around cold stars, it is also interesting to note that this planet is 
orbiting a relatively metal poor star, whereas I Santos et al.| ( [2003] ), 
Fischer & Valenti| ( |2005] ), and |Mayor et al.| ( |2011| ) have shown 
they have an even lower occurrence rate of gas giants for any 
orbital period. 

WASP-80b's equilibrium temperature will be around 800 K 
(for an albedo of 0.1). The planet- to- star contrast is favourable 
for future observation of the emission spectrum of the planet, 
because it is hosted by a star ~ 1500 K colder than the usual tar- 
gets. Furthermore, the near 3% depth of the transit makes this 
gas giant one of the most suitable targets for transmission spec- 
troscopy. 



Nota Bene We used the UTC time standard and Barycentric 
Julian Dates in our analysis. 
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